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Abstract. We report a study of microstructure and its quasi-periodicities of three pulsars at 1.65 GHz with 
the S2 VLBI system at a resolution of 62.5 ns, by far the highest for any such statistical study yet. For PSR 
B1929+10 we found in the average cross-correlation function (CCF) broad microstructure with a characteristic 
timescale of 95 ± 10 /is and confirmed microstructure with characteristic timescales between 100 and 450 ^s for 
PSRs B0950+08 and B1133+16. On a finer scale PSRs B0950+08, B1133+16 (component II) and B1929+10 
show narrow microstructure with a characteristic timescale in the CCFs of ~ 10 fis, the shortest found in the 
average CCF or autocorrelation function (ACF) for any pulsar, apart perhaps for the Crab pulsar. Histograms 
of microstructure widths are skewed heavily toward shorter timescales but display a sharp cutoff. The shortest 
micropulses have widths between 2 and 7 /^s. There is some indication that the timescales of the broad, narrow, 
and shortest micropulses are, at least partly, related to the widths of the components of the integrated profiles 
and the subpulse widths. If the shortest micropulses observed are indeed due to beaming then the ratio, 7, of the 
relativistic energy of the emitting particles to the rest energy is about 20 000, independent of the pulsar period. 
We predict an observable lower limit for the width of micropulses from these pulsars at 1.65 GHz of 0.5 /xs. If 
the short micropulses are instead interpreted as a radial modulation of the radiation pattern, then the associated 
emitting sources have dimensions of about 3 km in the observer's frame. For PSRs B0950+08 and B1133+16 
(both components) the micropulses had a residual dispersion delay over a 16 MHz frequency difference of ~ 2 fj,s 
when compared to that of average pulse profiles over a much larger relative and absolute frequency range. This 
residual delay is likely the result of propagation effects in the pulsar magnetosphere that contribute to limiting 
the width of micropulses. No nanopulses or unresolved pulse spikes were detected. Cross-power spectra of single 
pulses show a large range of complexity with single spectral features representing classic quasi-periodicities and 
broad and overlapping features with essentially no periodicities at all. Significant differences were found for the 
two components of PSR B1133+16 in every aspect of our statistical analysis of micropulses and their quasi- 
periodicities. Asymmetries in the magnetosphere and the hollow cone of emission above the polar cap of the 
neutron star may be responsible for these differences. 

Key words, pulsars: general - pulsars: radio emission, microstructure - methods: data analysis - methods: obser- 
vational 

1. Introduction recorded at different times are therefore related to differ- 
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ing neutron star. The radial modulation is likely related 
to plasma bunching and linked to the elementary emis- 
sion mechanism. In this model the spectrum of the radio 
emission is a function of the radial distance from the neu- 
tron star, and the beam width is frequency dependent. 
High frequency radiation is emitted closer to the neutron 
star and the beam is narrower, low frequency radiation is 
emitted further out and the beam is broader, reflecting 
the opening of the polar magnetic field lines. The study of 
pulsar intensity fluctuations has largely the goal of prob- 
ing on the one hand the geometrical characteristics of the 
pulsar emission beam and its underlying magnetospheric 
structure and on the other hand the elementary emission 
mechanism. 

Pulsar radio emission is known to exhibit fluctuations 
over a broad range of timescales. Av erage pulse pr ofiles 
can have up to seven components (Kramer 1994) and 
together with their frequency dependent widths reflect 
best the underlying geometrical structure of the magne- 
tosphere (Rankin 1983J). Every individual pulse is com- 
posed of one or several separate subpulses. In general, 
the subpulses fluctuate strongly within a single pulse and 
from pulse to pulse but have stationary characteristics and 
characteristic widths computed from their autocorrelation 
functions that are 96% correlated with the width of the 



strongest component of the average pulse profile (Bartcl 
et al. |1980| ). Like the average pulse profiles they most 
likely also reflect the geometrical structure of the mag- 
netosphere (Bartel et al. 198C). The subpulses in their 



turn are often composed of micropulses or microstructure 
with typical timescales of about hundred to a few hundred 
microseco nds, or several t enths of a degree in pu lsar longi- 
tude (e.g., Hankins 1972 ; Kardashev et al. 197S ). In a few 
cases still much faster but well resolved individual fluctu- 
ations were recorded, for inst ance with a timescale d own 
to 2. 5 fis for PSR B1133+16 flBartel fc Hankins||l982| , see 
also Bartcl 1978J ), the fastest fluctuations found for any 
pulsar apart from the Crab pulsar. For the latter pulsar 
sporadic giant pulses were observed which were stil l unre- 
solved at a time resolution of 10 ns ( Hankins 200C| ). 

For the broader micropulses, quasi-periodic structures 
were found in the very first studies of single pulses with 
sufficiently high time resolution. Hankins| ( |l97l| ) found 
many examples of regularly spaced micropulses with pe- 
riods of 300 to 700 /i s for P S R B 09 50+08 at a fre- 
quen cy of 1 1 1.5 MH z. |Backcr| fll973| ), [Boriakofl] ( [1976D 
and |Cordes| ( |l976a| ) have shown that PSR B2016+28 
has quasi-periodic microstructure with periods rang- 
ing from 0.6 to 1.1 ms at a frequency of 430 MHz. 



Soglasnov et al. (1981, 1983) analyzed the statistics of 
quasi-periodicit ies for PSRs B0809 +74 and Bl 133+16 
at 102.5 MHz. |Cordes et aL] ( |1990| ) studied five pulsars 
(including PSRs B0950+08 and B1133+16) with quasi- 
periodic microstructure at several radio frequencies. They 
concluded that there are no preferred periods for quasi- 
periodicities that are intrinsic to a given pulsar and that 
there is no frequency dependence of the micropulse width 
and the characteristic period of the quasi-periodicity. 



Lange et al. (1998) studied seven bright pulsars, in- 
cluding our three, at 1.41 and 4.85 GHz with a time reso- 
lution between 7 [is and 160 /is. They did not find notable 
differences of microstructure parameters at different fre- 
quencies. 

It is the topic of this paper to investigate micropulses 
and their quasi-periodicities and to help to understand 
whether they reflect the longitudinally modulated emis- 
sion pattern and the geometry of the magnetosphere (e.g., 
Bcnfordj 1977) or instead are more effected by the radially 
or intrinsic ally temporally modulated p ulsar emission pat- 
tern (e.g., |Hankins| [J972| ; |Cordesj |l98l|) . The shortest mi- 
cropulses observable in a pulsar in particular may harbor 
essential clues about the nature of the emission fluctua- 
tions in pulsars. 

To achieve a high time resolution one must digi- 
tally record the pulsar signal before detection with sub- 
sequent dispersion removal processing as originally de- 
scribed by Hankins ( 1971 ). Previous studies were based 
mainly on observations made with a time resolution of 
several microseconds to several tens of microseconds. A 



review of microstructure research is given by Hankins 
( 1996| ). In this paper we present a statistical analysis 
of the properties of microstructure for PSRs B0950+08, 
B1133+16, and B1929+10 at 1 650 MHz with a time reso- 
lution of 62.5 ns, the most extensive such analysis yet for 
any pulsar and with one of the highest time resolutions 
ever used. 



2. Observations 

The observations were made with the NASA Deep Space 
Network 70-m DSS43 radio telescope at Tidbinbilla, 
Australia. PSRs B0950+08 and B1133+16 were observed 
on 10 May 2000 and PSR B1929+10 on 24 April 1998. 

The data were recorded cont i nuously with the S2 
VLBI system (|Cannon et al.| |l997| |Wictfcldt ct~aL| |l998|) 
in the 2-bit sampling mode in the lower sideband from 
1 634 to 1 650 MHz and the upper sideband from 1 650 
to 1 666 MHz. Left circular polarization was recorded for 
both frequency channels. The observations were made in 
absentia which is more typical for VLBI observations. In 
general, pulsar observations with the S2 VLBI system can 
be made at any of the ~ 30 radio telescopes worldwide 
which are equipped with such a system, in the same way 
VLBI observations are made without the need for the in- 
vestigator's presence. In effect, a dedicated pulsar back- 
end at the observing station is replaced with a software 
package on the workstation at the investigator's home in- 
stitution. 



3. Data reduction 

The tapes were shipped to Toronto and played back 
through the S2 Tape-to-Computer Interface (S2-TCI) at 
the Space Geodynamics Laboratory (SGL) of CRESTcch 
on the campus of York University. The S2-TCI system 
transfers the baseband-sampled pulsar data to computer 
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Table 1. Pulsar characteristics and fixed parameters: P is the pulsar period; DM is the dispersion measure for a 
dispersion constant: ad = 2.4f x f CP 16 cm~ 3 pc • s; St ca x is the computed time delay between the two frequency 
channels based on DM, the standard errors are < O.f /is; N is the approximate total number of pulsar periods 
observed; N p is the number of selected strong pulses used to compute the CCFs. For PSR B1929+10 a distinction is 
made between the number used to compute the average CCF and the number (in parentheses) used to compute the 
individual CCFs for the analysis of individual pulses; N^ is the number of microstructure features revealed; N^-qp 
is the number of detected quasi-periodicities; AT no _ /x is the number of smooth structureless pulses; T is the number of 
samples for the duration of the pulse window used to calculate the CCF; SNR is the signal-to-noise ratio of the mean 
pulse intensity in the ON-pulse window averaged over all pulses selected for processing. It is a measure of the relative 
increase of the antenna temperature in the ON-pulse window. For the specific definition of SNR, see text. 
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files stored on hard disks and makes the data available 
for off-line analysis by a SUN Workstation computer. A 
similar use of the S2-TCI system for processing observa- 
tions of PSR J0437— 4715 was reported by Kcmpncr et al 

(S). 

The S2-TCI system enabled us to transfer the data 
stream to disk in a piecewise manner, the size or dura- 
tion of each piece of data stream being limited by the disk 
storage capacity. In general, we transferred a set of about 
10 minutes of data to disk at a time. In a first step we 
detected the recorded signal by squaring and averaging it 
with a time constant of about 100 /xs. Then we determined 
the phase of the ON-pulse window and selected strong 
pulses. For the selection we computed the signal-to-noise 
ratio, SNR = ifH , where (I on ) is the mean intensity in 
the ON-pulse window after subtraction of the mean in- 
tensity in the OFF-pulse window, (I s). The SNR so de- 
fined corresponds to the relative increase of the antenna 
temperature in the ON-pulse window and is therefore rel- 

_sma.ll even for strong pnlsps Wp nsoH SNR = 1 



— 1 to —3 is equal to the current root-mean-square (rms) 
deviation (+ler and —la, respectively) of the signal. In 
order to preserve this condition during observations, the 
S2 VLBI data acquisition system (S2-DAS) has an auto- 
matic gain control (AGC). For pulsar observations with 
the S2-DAS it is preferable to switch off the AGC and 
instead use the manual gain control with the gain fixed, 
or if left switched on, to choose a sufficiently long time 
constant for the AGC loop. Either option has the advan- 
tage of preventing the sampler from experiencing sudden 
gain discontinuities inside a pulse window. For our obser- 
vations the AGC was inadvertently left switched on, but 
fortunately the gain was found to be constant inside the 
selected pulse windows in the majority of cases. Therefore 
the threshold level could relatively easily be adjusted dur- 
ing the analysis after the observations to reflect the larger 
voltage variations in the ON-pulse window. We changed 
the decoding values through the data records from ±1 and 
±3 to real values that correspond to the current ±l a lev- 
els in accordance with the technique developed by Jcnct 



atively 

as a f.lrLsholH fa pulse sdertkm The a pprmH ma t P tntal & Anderson! (|l99g). We computed the new levels from 



number of pulses, N, observed and the number of selected 
pulses, N p , are listed for each pulsar in Table |[ This ap- 
proach reduced the amount of data by a factor of several 
hundred and enabled us to carry out the subsequent signal 
processing more efficiently 

Having determined the pulse windows and selected 
strong pulses we further processed the recorded (raw) sig- 
nal by decoding the signal amplitude sampled in two bits. 
Two-bit sampling of the amplitude of Gaussian random 
noise is widely used in VLBI observations. The decod- 
ing is generally done using four levels with intege r values 



the quasi-instantaneous rms values of subsequent portions 
of the data records, each 100 /is long, to approximately 
match the dispersion smearing time across the 16-MHz 
bandwidth. 

The next step in our data processing routine was 
the removal of the dispersion caused by the interstellar 
medium. The predetection dispersion removal technique 



equal to —3, — 1, +1, +3 ( [Thompson et al.||l988[) . These 
values reasonably represent the signal while the threshold 
level where the sampler switches from 1 to 3 and from 



itself (Hankins 1971) consists of a Fourier transform of 
the decoded signal followed first by amplitude corrections 
for the generally non-uniform receiver frequency bandpass 
and phase corrections for the dispersion delay, and then 
by an inverse Fourier transform back to the time domain. 
In particular, for the phase corrections of the dispersion 
delay, 5$(is), at the observing frequency, v, we used 
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8$(v) = (5$(^ + Al/) 



2tt DM 



-AS 1 



Av 



(1) 



where Av <C vq, DM is the dispersion measure, and 
«d is the dispersion constant (see caption of Ta b le 1) . 
Equation (|l|) was derived from Hankins fc Rickctt ( 1975 ) 
by conserving only quadratic and cubic terms in their 
equation (5). The next term would give a correction of only 
about 10 ns for the dispersion smearing over the 16-MHz 
bandwidth at 1650 MHz for a typical value of DM=3.0 
pc/cm 3 for our pulsars, which is only a sixth of our time 
resolution. Also, we did not need to apply a Doppler shift 
correction for DM due to the Earth's motion, since even 
for the largest possible velocity component of 30km/s such 
a correction would be limited to a smearing across the 16- 
MHz bandwidth of only about 40 ns for any of our pulsars 
which is smaller than our time resolution. The values of 
DM used for the three pulsars are listed in Table El After 
the amplitude and phase corrections we set to zero all 
negative frequency components in the spectrum and then 
computed the inverse Fourier transform. The sum square 
of the real and imaginary components of the complex time 
series r epresents the disper sion-free detected signal of the 
pulsar flBorn fc Wolj|l964[) . 

Finally, we calibrated the intensities. Their scale is 
based on the system temperature which was recorded in 
a log file. Noise fluctuations of the unsmoothed detected 
signal of la in the OFF-pulse window correspond to 50 Jy. 
Dynamically changing the conversion threshold levels dur- 
ing the ON-pulse and OFF-pulse windows permitted us to 
avoid a possible bias in the decoding of the two-bit sam- 
ples. 

4. Results 



4.1. Individual pulses 

In Figure [j] we display one example of a strong pulse for 
the upper (U) and lower (L) sideband for each of the three 
pulsars. The pulses are shown with a time resolution of 
32 fis. It can be seen that the pulses vary in intensity on 
a timescale down to almost the resolution interval. The 
intensity variations are strongly correlated between the 
upper and lower frequency bands. 

More rapid variations also appear to be present 
on first sight. For example, one very strong pulse of 
PSR B1133+16 has a mean flux density averaged over 
the pulse window of about 500 Jy, while the strongest 
isolated spike observed in the lower sideband has a flux 
density in excess of 6 000 Jy. A short 4-/zs portion of this 
pulse is shown in Figure Bwith the original time resolution 
of 62.5 ns. For such short intensity variations in time no 
correlation is apparent between the data of the two side- 
bands even if allowance is made for a possible variation 
of DM over a conceivable range of values. It is important 
to investigate whether these spikes and other strong un- 
resolved ones are related to the physical emitters or are 
just statistically insignificant noise fluctuations of much 
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Fig. 1. An example of a strong pulse at 1650 MHz (U) 
and 1 634 MHz (L) for each of the three pulsars. The dis- 
persion smearing was removed, but the dispersion delay 
between frequency channels was not removed. The inten- 
sity fluctuations were smoothed over a 32-/xs time interval. 

broader pulse structure that increased the ON-pulse an- 
tenna temperature and led to the strong variations. 

In the remainder of the paper we first discuss the 
properties of the microstructure with a width of ~ 1 to 
~ 500 /is. Then we present a search for structure on sub- 
microsecond timescales, and finally we report on a study 
of the statistics of microstructure quasi-periodicities of the 
micropulses of the three pulsars. 

4.2. Microstructure 

4.2.1. The typical microstructure width from the 
average cross-correlation function 

We determined the typical width of microstructure in a 
large number of single pulses from the width of the central 
spike of the averaged cross-correlation function (CCF). 
We chose CCFs for our analysis instead of the more 
widely used autocorrelation functions (ACFs) to avoid the 
noise spike at zero lag of the latter which could cause 
confusion in the detection of structure with the short- 
est timescales. More precisely, for each pulsar we selected 
N p (Table [jj) strong pulses, computed for each of them 
the CCF, i?i i2 (r), between the signals, Ii(t) and h(t), in 
the lower and upper sideband frequency channels, respec- 
tively, at lag r with 



#i,2(t) = [R hl (0)R 2 ,2(0)}- 1/2 J2h(.t)I 2 (t + T) 



(2) 



(=1 
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Table 2. Observed pulsar parameters: <5i bs is the observed time delay between the two frequency channels; <5i bs — St ca \ 
is the difference between the observed and computed time delay for the two frequency channels; T M _broad and T M „ narrow 
are the characteristic timescales of the microstructure determined from the average cross-correlation functions (see 
text); T^_ s hortest is the shortest width found for the observed micropulses ; ti/ 2 is the FWHM of the average profile or, 
in case of PSR B113 3+16 (I), of the com ponent in the pro file (taken from [Bartel ct al] ( |1980| )). For PSR B1133+16 (II) 
no width is given in Bartel et al, (198C); the width is from Kramer (1994) interpolated between his values of a FWHM 

(1980), is a measure of 



also from 



Bartel et al 



of a Gaussian at 1.41 and 4.75 GHz to 1.65 GHz. The parameter, 

the typical subpulse width determined as the 50% width (ignoring the zero-lag spike) of the average ACF of single 

pulses (for PSR B1133+16 r m = 2.9 ± 0.1 ms, no distinction is made between the two components) 
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Fig. 2. A short portion of a very strong pulse from 
PSR B1133+16 plotted with a 62.5 ns time resolution. 
The upper curve corresponds to the 1 650 MHz frequency 
channel, and the lower one corresponds to the 1 634 MHz 
frequency channel. The dispersion smearing and time de- 
lay between the frequency channels were removed. The 
original 2-bit quantization of the voltages is not visible 
because of the dispersion removal. 



Here t is the sample number at a particular time within 
the single pulse, T (see Table [j]) is the total number of 
samples chosen for the pulse window, and I\ (t) , 1% (t) rep- 
resent ON-pulse intensities with the OFF-pulse mean lev- 
els subtracted. 



We next averaged the CCFs for each pulsar and then 
corrected the averaged CCF for the receiver noise by di- 
viding it by the factor e (see Appendix IaI): 



1 



Coffi cr ff 2 



1 



l'°ni) . \-'on 2 / 
(ioffi) (J ff 2 ) 



(3) 



The quantities (loin), (-^on 2 ) represent the mean ON-pulse 
intensities inside the selected pulse window over N p pulses 
for each of the bands, 1 and 2, after subtraction of the 
equivalent mean OFF-pulse intensities, (/offi), (IoS 2 )- The 
quantities, a oni , cr on2 , a a 1 and <7 ff 2 are the rms devia- 
tions of Ion! , Ion 2 j Joffi , and I s 2 , respectively, again aver- 
aged over i V p pulses . This expression corresponds to equa- 
tion (20) in |Rickctt| ( |l975[) . 

Figure y shows the CCFs smoothed to a time-lag res- 
olution of 1 fis and shifted so that their time-lag origins 
correspond to the delays 5i ca i, computed from the dis- 
persion measure, DM (see Table []]). The left column 
shows the CCFs on a relatively large time-lag scale while 
the right column displays only the central part of the 
CCFs. The maxima of the CCFs are close to the value of 
0.5 pre dicted by the Amplitude-Modulated Noise (AMN) 
model ( Rickett 1975 ). The relatively small discrepancies 
may be due to small and possibly insignificant errors of 
the correction factors for the CCFs. 

For PSRs B0950+08 and B1929+10 a sharp spike can 
be clearly seen in the CCFs in the left column, indicating 
correlated microstructure in the two bands with a well de- 
fined characteristic timescale. For PSR Bl 133+16 a cen- 
tral spike is less pronounced in either of the two compo- 
nents at this frequency, but also existent. 

We determined the characteristic timescale of the typ- 
ical width of the microstructure from the intersection of 
the linear slope of the central narrow part (excluding the 
most inner section, see below) and the more slowly de- 
creasing broader part of the CCFs by linear least-squares 
fits. We list the values with the corresponding statistical 
standard errors as T„_broad in Table H. 
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Fig. 3. The average cross-correlation function (CCF) for PSR B0950+08, the first (I) and second (II) component of 
PSR B1133+16, and PSR B1929+10. The CCFs were calculated for the unsmoothcd ON-pulse intensities recorded in 
the conjugate 16-MHz bands and then smoothed over a time interval of 1 /is. The right column of the plots represents 
the very central portion of the CCFs shown in the left column. 
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Fig. 4. Examples of individual CCFs for PSR B1133+16 (component II) in order of increasing complexity. The left 
plot shows microstructure with one notable width only. The middle and the right plots show microstructure with two 
and three widths. The scale of the correlation coefficients is not corrected for receiver noise. 



Our value for r M _broad for PSR B0950+08 is com- ^^ 



patible with earlier measurements between 130 and 



obs 



200 [is (Pickett] [1975] |Cordcs fc Hankins| |1977| ; [Hankins 



& Bor%kofj|l978|; |Lange et al.||l998| ). For PSR B1133+16 
microstructure was reported in the range of ~ 340 to 



650 /is (Hankins 1972 ; Ferguson et al. 1976 ; Ferguson 



fc Seh]adakisj |l978j ; lCordcs| |l976a| ; popov et al.| |1987| ; 



Lange et al. 1995 ), comparable to our value of compo- 
nent I but three to sixfold larger than that for component 
II. For PSR B1929+10 the width of microstructure has 
never bee n measured in the a verage ACF or CCF before. 
However, Lange et al. ( 1998 ) reported on the detection 
of structure on timescales around 150 [is close to their 
effective time resolution of 70 /is. 

In the right column we display the central portions of 
the CCFs. For PSR B1133+16 (I) no additional abrupt 
sharpening of the CCFs can be seen. However for PSRs 
B0950+08, B1133+16 (II), and B1929+10 an additional 
central spike was found, indicating particularly short mi- 
crostructure with characteristic timescales, T M _ narrow , of 
order 10 /is (see Table 0). Such short microstructure has 
never been seen in the average ACF or CCF of any pulsar. 
However, for the giant pulses from the Crab pulsar, mi- 
crostructure was found in the ACF of a single pulse with 



the w idth at the 25% level of the ACF of < 1 /is (jHankins 
2000| ). 



4.2.2. Time delay of microstructure from the average 
cross-correlation function 

A close inspection of the inner part of the CCFs shows that 
the peak does not always occur exactly at the expected 
time delay. To quantify the discrepancy we measured the 
center of symmetry of the top part of the central portion of 
the CCFs in Figure y and list it as <5i bs in Table p[ To de- 
termine th e stan dard errors, A5i bs, we followed Chashe 
k Shishov| ( |1975[) and used 



4r M 
SNR 



(N P T) 



-1/2 



(4) 



where r^ is the microstructure timescale, r^-broad, listed 
in Table g, and SNR is the signal-to-noise ratio of the 
intensity of the average single pulse integrated over the 
pulse window (Table |l|) . 

For PSR B1929+10 the observed time delay is equal 
within A(5£ bs to the expected time delay. However, for 
PSR B0950+08 and the first and second component of 
PSR B1133+16 there are significant discrepancies. The 
differences, 6i bs — St ca .i, are about —2 /is (see Table J2|) 
in the sense that micropulses at the higher frequency ar- 
rive slightly later than expected in comparison to the mi- 
cropulses at the lower frequency. If interpreted as a disper- 
sion measure difference, they correspond to values of DM, 
1 to 2% smaller than those listed in Table UL However, 
as we will discuss in Section pi we do not think that the 
values of DM in Table [l] need to be corrected. 

4.2.3. The distribution of microstructure widths from 
individual cross-correlation functions 

In order to determine the distribution of microstructure 
widths for individual pulses, we computed the CCF for 
each of the N p selected strong single pulses. Note that for 
PSR B1929+10 the number in Table [j] is given in paren- 
theses, different from the number of pulses used to com- 
pute the average CCF. An interactive "TV" -based task 
was used to display the individual CCFs in a variety of 
timescales. With the help of a TV-cursor, the microstruc- 
ture width was determined, again, from that point of the 
CCF where the central steep slope flattens out. 

The majority of the selected strong pulses showed dis- 
tinct short-term microstructure in the individual CCFs. 
Three examples of individual CCFs are shown in Figure [| 
for PSR B1133+16 (component II). The left plot shows 
a CCF with an individual spike only, the middle and 
the right plots show more complex CCFs with two and 
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Fig. 5. Histograms of timescales of detected microstructure in the individual CCFs of single pulses for PSRs B0950+08, 
B1133+16 (components I, II), and B1929+10. The insets display histograms for the shortest timescales. 



three different widths for the microstructure. For PSRs 
B0950+08 and B1929+10 very often several, up to five, 
different structural scales were found in the same pulse. 
In these cases several values for the microstructure width 
were obtained from one pulse. The total number of values 
for the microstructure width is listed for each pulsar (and 
each component in case of PSR B1133+16) in Table H as 
Nfj,. A number of single pulses had smooth CCFs with- 
out any distinct microstructure features. That number is 
given as A no _ A , in Table pi. Only one from 225 pulses of 
PSR B0950+08, or less than 0.5% had no microstruc- 
ture. For PSR B1929+10 about 3% of single pulse had 
no microstructure, while for PSR B1133+16 about 17% 
and 38% of pulses had no microstructure in case of the 
first and the second component, respectively. This find- 
ing indicates even higher percentages of si ngle pulses with 



microstructure than repo rted earlier (e.g., Smirnova et al 



1994| ; [Langc et al] |1998| ) , if the results also hold for the 



weaker pulses which had to be ignored in our analysis be- 
cause of sensitivity reasons. 

The histograms of the microstructure widths are pre- 
sented in Figure g. For all three pulsars they are skewed 
towards shorter widths. In general, they show a moderate 
rise from broader to shorter width starting at ~ 200 /is 



for PSRs B0950+08 and B1929+10 and - 800 /is for 
the two components of PSR B1133+16. The rise becomes 
markedly sharper at ~ 40 /is for PSRs B0950+08 and 
B1929+10 and at ~ 200 /is for each of the components of 
PSR B1133+16, or at - 0.08 and - 0.02% of P, respec- 
tively. The histograms peak at ~ 20 to 30 /is for PSRs 
B0950+08 and B1929+10 and < 50 /is for the two com- 
ponents of PSR Bl 133+16. Towards shorter widths the 
histograms of PSRs B0950+08, B1133+16 (I, see inset of 
plot), and B1929+10 display a sharp cutoff at a width of 
- 10 /is, and of PSR B1133+16 (II) at < 10 /is (not vis- 
ible in plot). The shortest widths measured are 7 /is for 
PSR B0950+08, 6 /is and 2 /is for PSR B1133+16 (com- 
ponent I and II, respectively), and 5 /is for PSR B1929+10 
(Table g). The cutoff width and the shortest widths mea- 
sured are greater than a) the smoothing interval of 1 /is 
used in the analysis of the individual CCFs and b) the 
expected scattering time for either of the pulsars, and 
are therefore intrinsic properties of the pulsars. Again, 
there is a difference for the first and second component 
of PSR Bl 133+16. Micropulses with a width < 50 /is 
occur almost twice as frequent in component II than in 
component I. 
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4.3. Short-term noiselike intensity fluctuations 

Our average CCFs as well as our individual CCFs do not 
show any microstructure features with a sub-microsecond 
timcscale. On the other hand, as was mentioned in 
Section 4.1, strong pulses when plotted with the highest 
time resolution contain bright sub-microsecond fluctua- 
tions (see Figure ||). Also, Sallmcn ct al. (199E) reported 
to have found intensity fluctuations from the Crab pulsar 
which were still unresolved at their highest time resolu- 
tion of 10 ns. Is it possible that "nanopulses" exist in our 
data with a width not much larger than our highest time 
resolution of 62.5 ns? They may perhaps not be visible 
in our histograms with more than 100 times wider bins. 
They could perhaps also be largely uncorrelated for the 
two bands and therefore not apparent in the CCFs. To in- 
vestigate the significance of our fast intensity fluctuations 
we have to compare their statistics with the statistics of 
noise. 

We use two approaches in our data analysis: 1) com- 
putation of short-term ACFs with a time-lag resolution 
of 62.5 ns for ON-pulse and OFF- pulse windows, and 2) 
comparison of the distribution of the intensities of short- 
term (62.5 ns) fluctuations ON-pulse with the distribu- 
tion of such intensities OFF-pulse and also with the x 2 ~ 
distribution for thermal noise. 



submicrosecond micropulses, if present at all, do not con- 
tribute to the total power of pulsar signal variations by 
more than 0.3 %. 



4.3.2. Distribution of intensities 

We also compared the ON-pulse with the OFF-pulse in- 
tensity distributions and the ^-distribution for thermal 
noise for the data with the highest resolution of 62.5 ns. 
For each 8-/xs interval with T = 256 samples, the rms de- 
viation ct and the quantity, (I(t) — (I))/a, were computed. 
The results obtained for three extremely strong pulses of 
PSR B1133+16 (component I) are presented in Figure |[ 
Both distributions are fairly well fit by the theoretical 
curve for the ^-distribution. The deviations at large am- 
plitudes of both the ON-pulse and OFF-pulse distribu- 
tions from the theoretical line can be contributed to effects 
of the two-bit sampling. In general, no notable differences 
are apparent between the ON-pulse and OFF-pulse distri- 
butions. ON-pulse and OFF-pulse intensity fluctuations 
much shorter than 8 fj,s therefore have the statistics of 
thermal noise, consistent with the AMN model (Rickctt 
|1975D . The bright unresolved intensity spikes displayed in 
Figure || are therefore insignificant statistical outbursts. 
No nanopulses were found. 



4.3.1. Short-term ACFs 

Short-term ACFs, R(t), with 



*w = 7^ I>(*)-wro +?-)-</>] 



(5) 



were computed for single pulses from many short adjacent 
intervals of the data with detected signal, I(t), at sample 
number, t, both in ON-pulse and OFF-pulse windows. The 
short-time intervals were each 8 fis long and contained 
T = 256 samples. The instantaneous mean level (/) used 
for the subtraction is given as (I) = J2t=i ^W/^ 1 - For each 
single pulse m = 256 ACFs were computed for the ON- 
pulse window covering an interval of 2.048 ms centered 
close to the maximum of the selected pulses. The same 
number of ACFs was also computed for the OFF-pulse 
window of each single pulse. Note, that the definition of 
the ACF used here is not quite equivalent to the definition 
of the CCF in the equation (||). 

In order to obtain sufficient sensitivity, we selected 
only very strong single pulses of number k: 41 for 
PSR B0950+08, 24 for PSR Bl 133+16, and 30 for 
PSR B1929+10. For each of the pulsars no systematic 
difference was apparent between the ON-pulse and OFF- 
pulse average ACFs above the 3ctacf level. For an aver- 
age ACF of white noise: ctacf = /7 ! , , or ~ 0.001 in our 
cases. The corrections to the ON-pulse ACFs for receiver 
noise (see equation (g)) are only 5 — 10%. Therefore, since 
no deviations were found between ON-pulse and OFF- 
pulse ACFs above the level of 0.003, we conclude that 



5. Quasi-periodic microstructure 

5.1. Analysis method 

Micropulse periodicities can be revealed, e.g., by analyzing 
ACFs, CCFs, or power spectra. Since our observations 
were conducted in two adjacent frequency bands, we again 
use CCFs and in addition cross-power spectra. 

The cross-power spectrum is given as {Si^iffj,)^ with 



SM 



T/2 

£ 

=-T/2 



Ria(r)W(T)exp(-j2irf lt T) 



(6) 



where /^ is the frequency of the microstructure pe- 
riodicity, Ri^(j) is the CC F (see equation (Bh ), and 
W(t) is Hanning's window flBlackman fc Tukcy 1958|) 



to smooth 

2-kti 



variations of Si^if/i), with W{t) 



used 

~ [l + cos : 

Figure [7| shows an example of quasi-periodic mi- 
crostructure within a single pulse of PSR B1133+16 (com- 
ponent I) for the upper and lower bands, the correspond- 
ing CCF, and its cross-power spectrum. The effect of the 
quasi-periodicity of the microstructure is clearly seen in 
the CCF and the cross-power spectrum. 

In the following analysis we will use the cross-power 
spectra only. That has at least two advantages compared 
with the traditional use of ACFs and CCFs. First, in case 
of only one strong and isolated spectral line or feature 
as displayed in Figure |7| the cross-power spectrum pro- 
vides directly numerical values for the frequency, /„, and 
width, A/ M , of the detected feature that corresponds to 
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Fig. 7. An example of single pulse intensities with quasi-periodic microstructure for PSR B1133+16 (component 
I) recorded in the upper (U) and lower (L) sidebands (left) and the corresponding CCF (middle) and cross-power 
spectrum (right). For plotting purposes only, the intensities were smoothed with a time constant of 8 /xs, and the CCF 
was smoothed with a time constant of 1 /is and displayed with an arbitrary scale. For the cross-power spectrum we 
plot only the frequency range above 2 kHz, since the large power spectral density from features at lower frequencies 
would otherwise determine the scale completely. 



the quasi-periodicity. Second, in case of more than one iso- 
lated strong feature, such features can be relatively easily 
distinguished in the cross-power spectrum but not so in 
an ACF or a CCF. We selected a feature as "detected" 
if its power spectral density was larger than, or equal to, 
6 rms in the spectrum on both neighboring sides of the 
feature, or if the power spectral density exceeded a reason- 
able threshold. The last condition was necessary to detect 
strong spectral features in the very low frequency range 
where the local rms variation could be overestimated be- 
cause of the frequent complexity of the spectral features. 
In case of approximately symmetric features, the width, 
A/ M , was determined as the FWHM. In case of complex 
features we interpreted the feature as a blend of several 
narrow features and estimated the FWHM of the dom- 
inant unblended feature by measuring the half-width at 
half-maximum intensity from the apparently unblended 
side of the feature to the peak and then doubled that 
width. 

In Table we list the total number, -/V m _qp, of strong 
isolated features found in the cross-power spectra of N p 
single pulses for PSRs B0950+08, B1133+16 (I, II), and 
B1929+10. 



5.1.1. Types of microstructure periodicities 

In Figure |8j we show different types of cross-power spec- 
tra for several selected single pulses of PSR B1133+16. 
Similar types were also found for the other two pulsars. 
All spectra were normalized by their total power, i.e. the 
sum of all harmonics was set equal to unity. 

Visual inspection shows that all spectral features be- 
low 10 kHz are substantially broader than the applied fre- 
quency resolution of about 100 Hz. The spectra there- 
fore differ substantially from random noise and reflect the 
properties of the pulsar emission. 



The spectra may be classified in four general categories 
in order of increasing complexity. The first category com- 
prises spectra dominated by only one symmetrical feature 
(a-e) . The second category comprises spectra with two or 
more isolated symmetrical features with their frequencies 
being approximately multiple integers of the lowest fre- 
quency at which an isolated feature could be identified (f, 
g). These two categories correspond to the "classic" quasi- 
periodicities. However they constitute only about 5% of all 
analyzed spectra, about the same percentage for all three 
pulsars. The third category covers the spectra with several 
symmetrical isolated features located at more or less ran- 
dom frequencies (h, i). About 30-35% of all spectra fall 
under this category. Finally, the fourth category comprises 
spectra with isolated features at random frequencies as in 
the third category, but where each isolated feature is not 
symmetrical but instead structurally more complex (j, k). 
This is the most numerous category; it contains about two 
thirds of all spectra of all three pulsars. 



5.2. Histograms of the microstructure period 

In Figure O we show the histograms of the microstructure 
period, P M (Pf, = l/ffi), for PSRs B0950+08, B1133+16 
(I, II) and B1929+10. There are some general similari- 
ties between them. The histograms are skewed towards 
smaller microstructure periods. They have a relatively 
narrow peak at about 300 /is for PSRs B0950+08 and 
B1929+10 and a broader peak at ~ 800 /is for com- 
ponent I of PSR B1133+16 and ~ 400 /is for compo- 
nent II, two to fourfold larger than r^-broad- The his- 
tograms fall off sharply towards smaller periods. The 
widths of the histograms are relatively large. The bulk 
of microstructure periods (75%) falls in the range of 0.1- 
1.0 ms for PSR B0950+08, 0.2-3.0 ms for component 
I of PSR B1133+16, 0.1-2.0 ms for component II of 
PSR B1133+16, and 0.2-0.7 ms for PSR B1929+10. For 
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Fig. 6. Histograms of intensities at the highest resolu- 
tion of 62.5 ns recorded ON-pulse (rectangles) and OFF- 
pulse (stars) for three extremely strong single pulses of 
PSR Bl 133+16 (component I). The dotted line corre- 
sponds to the x 2 -distribution for thermal noise. The ordi- 
nate gives the number of occurrences as a fraction of the 
total number of intensity samples of 393 216. The abscissa 
gives the intensity in units of a. 
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Fig. 8. Four different types of cross-power spectra (see 
text) for selected single pulses of PSR B1133+16 for any 
of its two components. Each spectrum is normalized by the 
total power, i.e. by the sum of the power of all harmonics. 
All plots shown in the figure start from some relatively 
high frequency. Low-frequency features are much more 
powerful and would determine the scale of the plots such 
that the high-frequency features would become almost in- 
visible. The spectra were smoothed by using Hanning's 
window, reducing the originally high resolution only by 
~ 25%. The noise level in the plots is much smaller than 
the features displayed. The spectrum (e) is a reproduction 
of the spectrum in the right panel of Figure pi. 



PSRs B0950+08 and B1133+16 for which microstructure 
results were published earlier, these range values agree well 
with those obtained by others at other frequencies with a 
smaller sample size (Hankins 1971; Soglasnov ct al. 1981, 



1983| ; |Cordes et al.||l990| ; |Smirnova ct al.||l994| )' 

Again, for PSR B1133+16 there is a clear differ- 
ence in the histogram of microstructure periods. Quasi- 
periodicities with the shortest periods are relatively more 
numerous for the second component than for the first com- 
ponent. 

5.3. Histograms of the Q-factor of the microstructure 
periods 

The microstructure periods were determined from the 
peaks of the individual features in the cross-power spectra 
at frequency / M . The widths, A/ p , of the features var- 
ied widely. To determine the sharpness of the features 
or the "fuzziness" o f the m icrostructure periods we fol- 
lowed Cordcs ct al] ( 1990J ) and defined the Q-factor as 
Q = /^/A/ M . Figure |10| shows the histograms of the Q- 
factor for each of the three pulsars. A large majority of 



spectral features have relatively low values of Q with the 
most frequent value of about 1.5 to 2. In 75% of all cases 
Q < 5 for PSR B0950+08 and each of the components of 
PSR B1133+16, and Q < 4 for PSR B1929+10. However, 
some relatively high values with Q ~ 10 were also found, 
but only for PSR B0950+08 and the two components of 
PSRB1133+16. 

5.4. Average cross-power spectra 

In order to investigate the general aspects of the cross- 
power spectra we averaged the individual cross-power 
spectra for each of PSRs B0950+08 and B1929+10 and 
each component of PSR B1133+16. We then normalized 
each of the spectra so that they have equal power spectral 
densities at frequencies higher than about 50 kHz where 
they are largely flat and identical to the equivalent spectra 
obtained for the OFF-pulse signal. In Figure nl] we show 
these spectra for all three pulsars and compare them with 
the OFF-pulse spectra. At frequencies below about 20 kHz 
the ON-pulse power spectra show a large excess above the 
noise power spectrum. This excess reflects the presence of 
microstructure fluctuation power with timescales down to 
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Fig. 9. Histograms of the microstructure periods, P^, for the quasi-periodic features detected in the cross-power 
spectra. 



~ 50 /is. The value corresponds to the timescales of mi- 
crostructure revealed through the analysis of individual 
CCFs. The low frequency part can be approximated by 
the power law, S oc f~ a , with the exponent a ~ 2 for all 
three pulsars. In comparison, the exponent for the OFF- 
pulse spectrum is about 0.5. Such an increase of power 
toward the low frequencies is typical for low-noise ampli- 
fiers. The noise with these characteristics is called excess 
noise, low-frequency noise, or 1// noise. Usually its expo- 
nent is in the range from 0.8 to 1.5. Our rather low value 
of the exponent is probably due to the properties of the 
cross-power spectrum in our analysis. Clearly, the analysis 
of the average power spectrum is another way to detect 
the presence of microstructure and its quasi-periodicities. 

6. Discussion 

We have presented a statistical analysis of the properties 
of microstructure for PSR B0950+08, the two components 
of PSR B1133+16, and PSR B1929+10, the most exten- 
sive such analysis with one of the highest time resolutions 
(62.5 ns) ever made. 

Microstructure with a characteristic broad timescale 
between ~ 95 and ~ 450 fis in the ACF or CCF was con- 
firmed for PSRs B0950+08 and B1133+16 and for the first 



time measured for PSR B1929+10. Microstructure with a 
characteristic short timescale of only ~ 10 fxs was found for 
PSR B0950+08, the second component of PSR B1133+16, 
and for PSR B1929+10, the first time such a short char- 
acteristic timescale was detected in the average ACF or 
CCF from any pulsar. Clearly microstructure can have 
broad as well as short timescales for the same pulsar. This 
conclusion is further supported by the histograms of the 
microstructure timescales which generally showed a mod- 
erately rising part toward shorter timescales, correspond- 
ing to the broad micropulses, and a sharply rising part 
corresponding to short micropulses. 

Different characteristic timescales for micropulses were 
previously found by others, albeit for much broader 



widths, poglasnov et al.| ( |1981| ) observed PSRs B0809+74 
and Bl 133+16 at 102.5 MHz with a time resolution of 
10 /is and found that a long timescale of > 500 /us and a 
shorter one of < 500 /us are present in the ACFs of these 
pulsars. Later it was shown for the same pulsars that mi- 
cropulses with a long timescale are clearly correlated over 



a bro ad frequency range of 70-102.5 MHz ( Popov et al 



1987J) while short timescale micropulses are not correlated 



over a narrow frequency range of 1.6 MHz at a center fre- 
quency of 102.5 MHz (Smirnova et al. 1986). In contrast, 
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Fig. 10. Histograms of the microstructure period Q-factor. 



in the present study we have found for PSRs B1133+16 
(component II) and B1929+10 a clear correlation between 
the short timescale microstructure in the two adjacent 
bands separated by 16 MHz or 1% of the center frequency. 

Despite several hours of observations, no nanopulses 
were detected for any of the pulsars. The shortest 
micropulse was found for the second component of 
PSR B1133+16 and had a width of 2 /is. In particular, 
the histograms of the microstructure widths for all three 
pulsars showed a sharp cutoff at 5 to 10 /is. The rare oc- 
currence of micropulses with short widths of only a few mi- 
croseconds but with flux densities strong enough to be de- 
tected is consistent with earlier results for PSR Bl 133+16 



at nea rly the sam e frequency by Bartel ( 1978 ) and Bartel 
& Han rins ( 1982 ) who also found only one or at best a few 
micropulses with such a short width after several hours of 
observations. Either micropulses with widths shorter than 
about 1 /is do not exist at our observing frequency in the 
three pulsars we studied or they are so weak or occur so 
infrequent that they have not yet been observed. 



In contrast to our pulsars, [Hankins; ( |2000| ) reported 
to have found extremely strong fluctuations for the Crab 
pulsar which are still unresolved with a time resolution of 
10 ns. It remains to be seen whether the difference in the 
shortest timescale for microstructure for our three pulsars 



on the one hand and for the Crab pulsar on the other 
hand is caused by the Crab pulsar's young age, small pe- 
riod, or general uniqueness. Observations of other pulsars 
with a similarly high time resolution as ours are needed to 
investigate further this aspect of microstructure research. 

Detailed inspection of the central portion of the aver- 
age CCFs showed that the time delay of the micropulses 
between the two sidebands is for PSR B1929+10 within 
the error as expected, but is for PSRs B0950+08 and 
B1133+16 1 to 2% shorter than expected. The latter dis- 
crepancy has never been observed before. Earlier investi- 
gations of a possible discrepancy were based on observa- 
tions with a coarser timing accuracy and a larger relative 
and a bsolute frequenc y difference of the two channels than 
ours (Boriakofi 1983; Popov et al 



1999|) 



1987; sallmen et al 



The quasi-periodicities observed for all three pulsars 
were in most cases only relatively weakly pronounced with 
Q-factors generally between 1.5 and 2 and cross-power 
spectra showing rather complex features. Only in 5% of 
all cases did we detect classic quasi-periodicities with sin- 
gle spectral features and their harmonics. In fact, it ap- 
pears that the complexity of the spectra almost contin- 
uously covers all of our four different categories. In this 
sense there is no fundamental difference between spectra 
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Differences in microstructure parameters for the 
two components of PSR B1133+16 were first reported 
by Cordes & Hankins (1977) in their study of polariza- 
tion properties of microstructure at 430 MHz with a time 
resolution of 8 /us. Also, Smirnova et al. (1994) found that 
in the range from 200 to 1 500 (is micropulses of compo- 
nent II more frequently had smaller widths than those of 
component I in their observations at 1.4 GHz with a time 
resolution of 150 (j,s. These results are at least qualitatively 
consistent with our result. 

Are the micropulses we observe due to a longitudi- 
nal emission modulation and therefore to a sweep of the 
beam past the observer, or due to a radial emission mod- 
ulation and therefore to the intrinsic time variability of 
the emission, or perhaps due to a combination of both? 
The average pulse profile clearly reflects the beam pattern; 
the equivalent profile width is proportional to P ( Lyne & 
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Fig. 11. Average cross-power spectra of the intensities 
at two conjugate 16-MHz bands for ON-pulse and OFF- 
pulse emission. The spectra were normalized to the largely 
flat part of the OFF-pulse spectrum between 50 kHz and 
1 000 kHz. 



of classic quasi-periodicities and spectra with no clear pe- 
riodicities at all. 

In this context the interpretation of microstructure 
quasi-periodicities in terms of modes of vibrations of neu 
tron stars (Boriakofi 1976; Van Horn 



1980; Hansen & 



Cioffi t_98C ) is not attractive although some modes fall 
well into the range of the observed quasi-periodicities (see 
McDcrmott et al. 1988). However, our values of the Q- 
factor are very small in comparison to about 1 000 ex- 
pected from stellar vibrations. Further, the essential con- 
tinuum of the degree of complexity of the cross-power 
spectra and the lack of any resonance spike in the his- 
tograms for the microstructure period, P^, make it very 
unlikely that the microstructure quasi-periodicities are 
caused by vibrations of the neutron star . Similar conclu- 
sions were drawn by |Cordes et al. (199C) as a result of a 
study of five pulsars at several frequencies. 

We find it most remarkable that the micropulses in 
the two components of PSR B1133+16 display different 
characteristics in every aspect we analyzed. The aver- 
age CCFs, the histograms of microstructure timescales, 
the histograms of microstructure period, and the average 
cross-power spectra are all significantly different for the 
two components. Most striking is the more dominant oc- 
currence of short micropulses in the second component 
which is reflected in the sharp spike of the central portion 
of its average CCF and the twofold higher frequency of 
occurrence in the histogram for < 100 /is widths. 



Manchester 1988 ). Subpulses in single pulses have a typ- 
ical width oc P 08 indicating that the millisecond struc- 
ture like the average pulse profile can largely be inter- 
preted as being due to the sweep of the beam past the 
observer (Bartel et al. 198C). Broad micropulses of sev- 
eral 100 [is duration were reported to have a somewhat 
weaker dependence on period (oc P - 5 ; Taylor et al. 1975; 



see also Ferguson 1977) reflecting perhaps already part 
of the elementary emission mechanism with its temporal 
variations. 

Three pulsars are not enough to test this relation since 
the scatter about the regression line is relatively large. 
However, we can compare our broad, narrow, and short- 
est microstructure widths with the component widths of 
the average profile, £1/2, and with the subpulse widths, r m , 
listed for the three pulsars in Table Q The ratios for the 
comparison are given in Table [3[ While the microstruc- 
ture widths relative to P can indeed vary within their cat- 
egories by several times their errors , they are constant 
within about twice the combined errors relative to £1/2 
and r m (with the exception of T M _broad relative to ti/2 for 
component I of PSR B1133+16). 

In other words, with the exception of the broad mi- 
cropulse widths for PSR Bl 133+16 (I), each of the differ- 
ent types of micropulse width scales within ~ 2er with the 
subpulse width and the average pulse component width. 
That may be an intriguing result. It therefore appears pos- 
sible that even the fastest observable intensity fluctuations 
still have a strong component of the longitudinal emission 
modulation and are at least partly due to the sweep of the 
beam. If so, we can compute the factor 7 as the energy rel- 
ative to the rest energy of the particles moving along the 
curved magnetic field lines and beaming in the direction 
of motion. For point sources, 



1 



7 



sin £ 27r sin £ 



-narrow 



(see Lange et al. 1998 ). Here, (ft is r ( 



(7) 



in radians 



relative to P, and C is the angle between the rotation axis 
of the pulsar and the line of sight to the observer. For the 
conservative assumption of £ = 90° the values of 7 differ 
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Table 3. Comparison of pulsar parameters and 7 factors. The parameters for comparison in columns 2-4 are taken 
from Table g. For PSR B0950±08 3% errors are assumed for ti/2 and r m . The parameter 7y_~ narr o W is derived for an 



angle, £, between the line of sight and the rotation axis of the pulsar of 90° and 7/i-narrow for an angle estimated from 
the polarization angle sweep across the bea m for the hollow-cone mod el, £ = 10.1, 55.0, 10.0 for PSRs B0950±08, 
B1133±16 (II), and B1929±10, respectively (Lyne & Manchester 1988). For the first two pulsars the values for ( are 
similar to those given by Rankin (1993). For PSR B1929±10 Rankin (1993) interpreted the viewing geometry very 
differently in comparison to Lyne & Manchester (1988), but did not give a value of £. The errors in columns 5-6 are 
"symmetrized" errors based on the errors of T^_ narrow only. 



PSRB 


7"/i — broad 


7~fi — narrow 


7~l-i — shortest 


(C=90°) 

1 (J, — narrow 


OV 4 — narrow 


C ' Tfi — narrow 




(HT 4 P ) 


(10" 5 P ) 


(10- 3 r m ) 






(km) 




(10- 2 ti /a ) 


(HT 3 t 1/2 ) 












(10~ 2 r m ) 


(10- 3 r m ) 










0950±08 


5.3 + 0.2 
1.5 ±0.1 
2.7 ±0.1 


5.5 ±1.2 
1.5 ±0.3 
2.8 ±0.6 


1.40 ± 0.04 


2 900 + 700 


16 400 + 4 000 


4.2 + 0.9 


1133+16 (I) 


3.6 ±0.3 
7.2 ±0.9 


: 




~ 




" 


1133+16 (II) 


0.9 ±0.2 
1.0 ±0.2 


0.9 ±0.3 
1.0 ±0.3 


- 


16 200 ± 5 000 


19 700 + 6 000 


3.3 + 0.9 


1929+10 


4.2 ±0.4 
1.6 ±0.2 
2.8 ±0.3 


4.0 ±1.3 
1.5 ±0.5 
2.7 ±0.9 


1.47 + 0.04 


4 000 + 1500 


23 000 + 9 000 


2.7 + 0.9 



by up to a factor of 5. However, if the angle ( is taken 
from estimates obtained from the sweep of the polarization 
angle across the pulse profile (Lyne & Manchester 1988| ), 



the values for 7 are much more uniform with a mean of 
about 20 000 and a spread of not more than a factor 1.4 
which is well within the errors. Perhaps such a high value 
for 7 is representative for the relativistic energy of the 
particles. 

Any longitudinal emission pattern of the beam on this 
small angular scale is likely related to the radial pattern 
through the radius to frequency relation and the curved 
open magnetic field lines of the polar regions. For instance, 
for PSR B0950+08 the subpulse ACF width i ncreases from 
4.2 m s at 2.7 GHz to 5.0 ms at 1.7 GHz ( |B artel ct aT 
1980]), reflecting the above frequency dependent geome- 



try. Since the subpulses can be approximated to first or- 
der as Gaussiansn, the 50% ACF width has to be scaled 
upwards by a factor of 1.4 to yield the FWHM of the sub- 
pulses. If we had observed the subpulses over a continuous 
bandwidth from 2.7 to 1.7 GHz, they would be smeared in 
time at the wings by about a third of the difference of the 
scaled widths or 0.37 ms. On average, the smearing time 



A detailed analysis showed that the average ACF of sub- 
pulses of a representative sample of pulsars, and by implication 
the average subpulse itself, is sharper at the top than either 
a Gaussian or a Lorenzian but flatter than a Gaussian and 



similar to a Lorenzian in the wings (Bartel et al. 1980) 



for a subpulse would be about 0.19 ms. Correspondingly, 
we could expect our micropulses observed over a total of 
the two bands of 32 MHz to be smeared in time on aver- 
age by about 6.1 us. For PSRs B1133+16 and B1929+10 
the expected smearing times based on the differences of 
the scaled ACF widths are 4.5 (see caption of Table S for 
PSR B1133+16) and 3.7 /is, respectively. The expected 
smearing times for the first and third pulsar are in good 
agreement with their shortest widths. For PSR B1133+16 
the ACF width refers to a weighted mean of the ACF 
widths of both components. The mean of our shortest 
widths for both components is 4.5 (is which equally well 
agrees with the expected smearing time for both compo- 
nents. 



In this context the residual delay of ~ 2 /is that we 
observed for the microstructure of the two bands for two 
pulsars can be interpreted as being due to the curved ge- 
ometry of the magnetic field lines and to the radius to 
frequency mapping. One might expect the delay to be a 
function of the phase of the micropulses in the pulse win- 
dow. Micropulses in the leading part of the pulse window 
would first appear at the lower frequency, micropulses in 
the trailing part first at the higher frequency. In case of 
PSR BI 133+16 for which we distinguished between a lead- 
ing and a trailing component we did indeed find that the 
lower frequency micropulses arrived first in the leading 
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component. However, we did not observe a sign reversal 
of the delay for the micropulses in the second component. 

This complication in the interpretation notwithstand- 
ing, it is conceivable that the smallest observable angular 
scale of the fluctuating radiation pattern of the beam is 
limited by the geometry of the opening cone of emission 
and therefore likely by the opening polar magnetic field 
lines, the particular radius to frequency mapping, and the 
bandwidth of the receiver. Shorter pulses in these pulsars 
may become observable by using a smaller bandwidth. The 
shortest pulses observable at 1.7 GHz would be limited by 
the rise time of the optimal filter to about 0.5 fin. Other 
pulsars may have a more favorable, weaker dependence of 
the ACF width on frequency and could show less angular 
smearing. A detailed study of the correlation and the pulse 
time of arrivals of micropulses at several close frequencies 
with a similar resolution of 62.5 ns could provide impor- 
tant clues about the geometry of the magnetic field lines, 
the elementary pulse structure, and in turn the emission 
mechanism. 

However, the alternate interpretation of the short mi- 
cropulses being due to a radial modulation of the radiation 



tions are envelope solitons or isolated wave packets and pe- 
riodic wave trains. The first solution may describe strong 
isolated micropulses while the second could describe the 
phenomenon of quasi-periodicity. 

Recently Weatherall ( 1998 ) conducted a numerical 
simulation based on nonlinear wave dynamics which pro- 
vides detailed solutions for the temporal behavior of pul- 
sar radio emission. The model predicts an intrinsic pulse 
width in the range of 1-10 ns (nanostructure). The charac- 
teristic timescale is expected to be frequency dependent, 
becoming longer at lower frequencies. 

Our analysis of pulsar radio emission with the res- 
olution of 62.5 ns has not detected any nanopulses. 
Nanostructure may in general be present in radio emis- 
sion as "shot noise" as suggested by Cordcs (1976b). 
Perhaps nanopulses might be observed as an infrequent 
phenomenon like the "giant" pulses from the Crab pulsar. 

However, no matter how short elementary pulses could 
intrinsically be, they are not only likely affected by angu- 
lar smearing as suggested for our pulsars above but could 
perhaps also be significantly affected while propagating 
through the magnetosphere toward the observer. Lyutikov 



pattern is not necessarily disfavored. In the context of this fc Parikh (200C) have studied the scattering and diffrac- 



interpretation the values for 7 could be easily two to three 
orders of magnitude lower and the typical radial length of 
the modulation ~ c • T^-narrow (see Table y). The sources 
of such short micropulses have brightness temperatures of 
~ 10 34 K. With the influence of beaming being excluded, 
the high degree of coherence of the modulation over the 
radial length would primarily account for the high bright- 
ness temperature. 

Several possible emission mechanisms have been pro- 
posed to explain radio emission of pulsars. In general, they 
can be classified into three groups: emissi on by bunches , 
plasma instabilities, and maser emission ( Mclroscj 1992J) . 
Theoretical predictions on the possible generation of ultra 
short timescale intensity fluctuations in the pulsar radio 
emission are based on nonlinear temporal models describ- 
ing the interaction of a high energy beam of particles with 



plasma wave packets in the pulsar magnetosphere. Asseo 



et al. ( 1990 ) suggested that a two-stream instability may 
result in the growth of strong plasma turbulences which 
would lead to a self-modulation instability, the creation 
of localized wave packets, and the generation of stable 
Langmuir soliton-like solutions. The latter would result in 
the so-called "Langmuir microstructures" with timescales 
in the range of several microseconds. Our shortest mi- 
cropulses have the predicted typical durations and could 
be interpreted in this way. 

Other interpretations of microstructure are based on 
nonlinear models of self-modulational instabilities of an 
electromagnetic wave with a large amplitude propagat- 
ing in an electron-positron plasma ( Chian fc Kcnne| 1983 ; 
Onishchcnkcl |l990| ; |Chian||l992t pangadhara et al.||1993|) 



The models demonstrate that the modulation instability 
can evolve to a nonlinear stationary state that results from 
the balance between nonlinearity and dispersion. The pos- 
sible stationary solutions of the respective nonlinear equa- 



tion of radio waves in the pulsar magnetosphere and found 
that the scattering delay can be as large as 10 /xs, easily 
explaining our cutoff widths in the histograms. The dis- 
persion delay of micropulses of about 2 fis we found for 
PSR B0950+08 and both components of PSR B1133+16 
could also be explained in their model as a propagation 
effect. 

The difference in the characteristics of microstructure 
for the two components of PSR B1133+16 remains puz- 
zling. If the two components define the longitudes where 
the hollow cone of emission is intersected by the line of 
sight, then the emission mechanism produces micropulses 
with characteristics as a function of azimuth about the 
axis of the hollow cone. Perhaps asymmetries in the po- 
lar cap magnetic field structure could cause azimuthally 
dependent plasma bunching properties and propagation 
effects that would result in the differences we observe. 

7. Conclusions 

The main results and the conclusions of our study can be 
summarized as follows: 

1. The S2 VLBI system provided a conceptually clearly 
defined continuous flow of data with a time resolution 
of 62.5 ns that could be transferred via the S2-TCI 
to small digital tapes for the analysis of pulsar radio 
signals. 

2. Broad microstructure was found for the first time in 
the average CCF (or ACF) of PSR B1929+10 and con- 
firmed for PSRs B0950+08 and B1133+16. 

3. Narrow microstructure with a typical width of only 
about 10 /is was discovered in the average CCF 
of PSRs B0950+08, B1133+16 (component II) and 
B1929+10. 
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4. Histograms of microstructure timescales show a steep 
increase toward shorter timescales followed by a sharp 
cutoff at about 5-10 /is. The shortest micropulse de- 
tected had a width of 2 /j,s. 

5. No unresolved nanopulses or pulse structure with a 
sub- microsecond timescale were found. This conclusion 
follows from the comparison of histograms of ampli- 
tudes for ON-pulse and OFF-pulse windows; it is con- 
firmed by the completely identical shapes of the high- 
frequency portions of the average cross-power spectra 
for ON-pulse and OFF-pulse emission; and it is sup- 
ported by the short-term ACFs which show that the 
pulsar signal is indistinguishable from pure receiver 
noise at timescales shorter than 8 /its. 

6. For PSR B1929+10 the dispersion delay for mi- 
cropulses over 16 MHz or ~ 1% of the center frequency 
is the same within one standard error as that for the 
average pulse over a much larger absolute and relative 
frequency spread. For PSRs B0950+08 and B1133+16 
the residual delay is about 2 /is with the micropulses 
at the lower frequency arriving earlier than expected. 

7. The cross-power spectra for single pulses display fea- 
tures with an almost continuous variety of complexity 
ranging from single features that correspond to the 
classic quasi-periodicities to multiple and broad over- 
lapping features with random positions in the spectra 
that correspond to essentially no quasi-periodicities at 
all. 

8. The power spectral density of the average cross-power 
spectra of single pulses decreases steeply toward higher 
frequencies of the microstructure quasi-periodicities 
and merges with the power spectral density of the 
OFF-pulse emission at a frequency that corresponds 
to the shortest microstructure quasi-periodicities. 

9. The statistics of the micropulses and their quasi- 
periodicities differ significantly for the two components 
of PSRB1133+16. 

10. It is conceivable from our data that the shortest mi- 
cropulses represent, at least in part, the shortest an- 
gular scale of the longitudinal radiation pattern of the 
sweeping beam with 7 ~ 20 000 for the radiating par- 
ticles. The particular combination of the geometry of 
the opening, curved magnetic filed lines, the radius to 
frequency mapping, and the bandwidth over which the 
radiation is observed limited the observable width of 
micropulses in our observations to ~ 5 fis. The short- 
est micropulses of our three pulsars at 1.7 GHz are 
predicted to have a width of ~ 0.5 /is. If the short mi- 
cropulses are instead due to a radial modulation of the 
radiation pattern the typical length of the pattern is 
~ 3 km. 
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Appendix A: Correction of the average CCF for 
the effects of receiver noise 

The effects of receiver noise on the average ACF of a pulsar 
signal were investigated by Rickett ( 1975| ). The ACF, r(r), 
of the pulse intensity in terms of observable quantities is 
given as 



(r(r)) = <r on (r))-<r off (r))(l + 2^) , 



where r on (r) and r e(r) are the ACFs for the ON-pulse 
and OFF-pulse windows each averaged over all pulses. 
Further, (/on) is the mean intensity in the respective win- 
dow again averaged over all pulses, but with (/ fif) sub- 
tracted. The normalized ACF, R(t), is then given as 



(R(r)) = 



(r(r)> 



(ron(r)) 



< r (°)> (r on (0))-(r off (0))(l + 2^ 



(foff( T )) being essential zero for t > 1/B, where B is the 
receiver bandwidth. 

This equation can be rewritten in terms of (i? on (T)) = 

(ron(T)) 



as: 



(R(r)) = 



(Ron(r)) 



where e is a correction factor: 



1 _ °" off I 1 1 oifon) 



with cr on = y/(r on (0)) and er off = y/(r on {0)). 

In this paper we used the CCF between adjacent 16- 
MHz frequency channels. It can be seen that for the CCF 
the correction factor e should be modified to: 



€= 1 



Coffi Ooff 2 



V*oni ) 1 (-'0112 ) 
(/offi) (/off 2 ) 



In fact, in our observations, the values of a on , a g, 
(/ on ), and (/off) in one channel never differed from the 
equivalent ones by more than a few percent in the other 
channel because of the proximity of the two frequency 
channels and the AGC operation. 
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